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I.  INTRODUCTION 

Recently  there  has  been  much  interest  in  the  properties 
of  quasiperiodically  driven  nonlinear  oscillators  of  vari¬ 
ous  types.  In  addition  to  the  more  well  known  dynamical 
behaviors  which  result  in  two-frequency  quasiperiodic, 
three-frequency  quasiperiodic,  and  chaotic  attractors, 
such  oscillators  have  been  shown  to  exhibit  a  new  type  of 
dynamical  behavior  leading  to  what  are  termed  strange, 
nonchaotic  attractors  (SNCA’s).  This  behavior  was 
demonstrated  analytically  and  shown  to  be  structurally 
stable  [1],  and  subsequently  verified  in  numerical  itera¬ 
tions  of  maps,  representing  dissipative  nonlinear  systems 
with  periodic  or  bistable  potentials,  driven  by  two  period¬ 
ic  forces  with  incommensurate  frequencies,  that  is,  quasi¬ 
periodic  forcing  [1-6].  Moreover,  a  SNCA  has  been  ob¬ 
served  also  in  the  map  representation  of  an  infinitely 
damped  system  [2,3]  which  shows  no  chaos  or  three- 
frequency  quasiperiodicity.  Though  some  early  doubts 
existed  [2],  the  SNCA’s  observed  in  these  systems  were 
later  shown  to  exist  on  a  subset  of  parameter  space  of 
measure  greater  than  zero  [2-5].  Since  the  SNCA’s  were 
numerically  demonstrated  in  the  resisti vely  shunted 
Joscphson  model  [2-4]  and  in  an  oscillator  with  the  stan¬ 
dard  quartic  potential  (Duffing  oscillator),  which  serves 
as  a  generic  representation  of  a  variety  of  realizable  bist¬ 
able  nonlinear  systems,  and  since  they  exist  on  a 
positive-measure  subset  of  parameter  space,  it  seemed 
reasonable  that  they  could  be  discovered  in  an  actual  ex¬ 
periment.  Indeed,  the  experimental  demonstration  of  the 


existence  of  a  SNCA  in  a  quasiperiodically  forced,  free¬ 
standing  beam  of  magnetic  material  whose  mechanical 
properties  can  be  externally  controlled  has  been  recently 
accomplished  [7,8].  In  addition,  extensive  numerical 
studies  on  a  model  of  the  magnetic  beam  apparatus  based 
on  the  parametrically  driven  Duffing  oscillator  have  been 
carried  out  [8]. 

The  emphasis  of  much  of  the  numerical  work  has  been 
on  the  development  of  criteria  by  which  SNCA’s  could 
unambiguously  be  characterized  and  on  studies  of  the 
route  to  chaos  in  two-frequency  quasiperiodically  forced 
systems.  SNCA’s  show  a  largest  Lyapunov  exponent 
which  is  nonpositive,  hence  the  term  nonchaotic.  Beyond 
that,  the  two  most  useful  characterizations  which  have 
emerged  are  based  on  the  power  spectrum  and  on  the  in¬ 
formation  dimension  Z>t  extracted  from  a  Poincare  sec¬ 
tion.  The  Poincare  section  itself  also  displays  differences 
from  those  of  chaotic,  quasiperiodic,  or  periodic  attrac¬ 
tors:  in  theory,  its  structure  is  everywhere  single  valued 
but  everywhere  discontinuous,  hence  the  term  strange.  In 
a  typical  system,  for  example,  the  quasiperiodically  drive  A 
pendulum,  the  two  Lyapunov  exponents  (of  the  reduced 
map,  created  by  sampling  the  two  variables  in  phase  with 
one  of  the  driving  functions)  are  found  to  be  zero  and 
smaller  than  zero,  respectively  [4].  Following  the 
Kaplan-Yorke  conjucture  [9],  the  information  dimension 
of  such  an  attractor  should  therefore  be  unity  [6],  The 
one  experiment  which  has  been  accomplished  [7]  mea¬ 
sured  a  somewhat  larger  value  /),  >1.3.  t  he  numerical 
expenmeiits  of  Ding,  Grebogi,  and  Oil  jo]  yielded  values 


of  £>,  sa  1.2  for  the  reduced  map  of  the  pendulum  and 
Z>,£50.9  for  a  quasiperiodically  forced  circle  map.  We 
show  here  in  Sec.  Ill  that  careful  measurements 
(  sa  100000  points  in  the  Poincare  section)  on  an  analog  > 
simulator  result  in  a  value  in  quite  good  agreement  with  : 
these.  £>,«  1.2,  which  can  be  compared  to  our  value  of 
£),  3  1.6  in  a  very  close  region  of  parameter  space  which  j 
is  chaotic.  The  following  additional  characterization  has  j 
been  developed.  From  the  power  spectrum,  the  number ; 
of  peaks  .V(a )  exceeding  a  threshold  amplitude  a  obeys  a  j 
simple  power-law  scaling  unlike  those  for  two-frequency  j 
quasiperiodicity  or  for  chaos  [2,3].  Specifically,  : 

withl<a<2  (la)  j 


for  the  SN'CA, 


FIG.  1.  A  Plot  of  Eq.  (4)  with  0—2,  which  was  the  potential  ’ 
used  in  this  work. 
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for  two-frequency  quasiperiodicity,  and 
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for  three-frequency  quasiperiodicity.  All  of  these  charac¬ 
terizations  have  now  been  investigated  numerically  in  the  : 
actual  experiment  and  in  the  analog  simulation  described  i 
here.  Some  of  them  have  been  used  in  numerical  studies 
on  SNCA’s  in  the  Van  der  Pol  and  other  oscillators 
[10-13].  Studies  of  the  quasiperiodically  forced  Duffing 
oscillator  were  accomplished  by  Wiggins  [14]  and,  in  an  j 
interesting 'experiment  with  a  mechanical  oscillator  : 
which  can  be  represented  by  the  generic  Duffing  system  : 
[14],  by  Moon  and  Holmes  [15].  In  neither  Ref.  [14]  nor  : 
Ref.  [15],  however,  was  strange  nonchaotic  dynamics  re-  >■ 
ported. 

In  this  paper,  we  investigate  the  existence  of  SNCA’s  : 
and  the  transition  to  chaos  both  by  digital  numerical  I 
simulation  and  by  measurements  made  on  an  analog 
simulator  of  a  version  of  the  pendulum  equation  which 
has  been  used  to  represent  the  radio-frequency-drivcn  su¬ 
perconducting  quantum  interference  device  (SQUID) 
with  inertia  and  finite  damping: 


x  +  k.x  ~  — 


dU(x) 

dx 


+  <7,sin<j,t  +  q-isifua1t  » 


(2) 


where  the  ratio  of  the  two  frequencies  is  the  golden  mean 
l+t/5 


(ji  2 

The  potential  used  in  this  work  was 


0. 

2rr 


cos2n-x  . 


(3) 


(4) 


and  for  all  [he  results  reported  here  we  took  0  -  2,  which 
results  in  the  tristablc  potential  shown  m  Fig.  |.  This 
•  entrusts  with  the  single  previous  experiments  [7,8] 
•vhe;  h  used  a  more  or  less  symmetric  bistable  potential, 
J?1 1,{  the  previous  numerical  work  which  used  either 
,  r  histablf  [H.rrith  ik  Our  moitvatift*  for  tins 


work  was  to  demonstrate  the  existence  of  SNCA’s  in  this 
particular  potential  by  experimental  measurements  on  a 
real  physical  system,  the  analog  simulator,  and  thereby  to 
demonstrate  the  practicality  of  an  experiment  on  a  sys¬ 
tem  with  a  similar  potential,  perhaps  a  r.f  SQUID.  Our 
results  further  illustrate  the  robust  properties  of  this  dy¬ 
namics  and  the  case  with  which  it  can  be  observed  in  real 
physical  systems. 

In  Sec.  II,  wc  describe  the  simulator  together  with  as¬ 
sociated  apparatus  and  the  method  of  obtaining  experi¬ 
mental  data  from  it.  The  results  from  the  analog  Simula- . 
tor  are  backed  up  by  numerical  experiments,  the  methods  : 
for  which  are  also  described  in  this  section.  In  Sec.  Ill, 
we  display  characteristic  results,  in  particular  those  for . 
the  information  dimension,  Lyapunov  exponent,  and  in-  j 
variant  density  for  both  chaotic  and  strange  nonchaotic 
attractors.  Measurements  of  representative  power  spec-  ( 
tra  are  displayed  and  analyzed  according  to  Eqs.  (1). 
Moreover,  we  demonstrate  that  a  SNCA  can  destroy  and 
replace  a  chaotic  attractor  under  the  action  of  a  single  | 
control  parameter.  We  display  this  result  in  the  style  of  a  ■ 
phase  transition,  using  the  Lyapunov  exponent  as  the  or¬ 
der  parameter.  In  Sec.  IV,  we  summarize  our  results  and 
suggest  future  experiments. 

II.  SIMULATION  TECHNIQUES 

Analog  simulators  for  various  applications  have  been 
described  before  [15]  in  both  antique  [16]  and  modern 
[17]  versions.  Our  simulator  of  the  dynamics  of  Eqs.  (2) 
and  (4)  is  shown  in  Fig.  2.  It  operates  over  a  dynamic 
range  of  ±  10  V  with  a  resolution  (nearly  an  order  of  mag¬ 
nitude  above  the  inherent  circuit  noise)  of  s  10  mV.  The 
voltages  x(t)  and  i(f)  were  digitized  in  pairs  at  12-bit 
resolution  by  a  Data  Translation  Model  DT2828  analog* 
to-digital  converter  (ADC),  which  had  a  throughput  of 
130  kHz,  and  input  to  a  Model  AT  personal  computer 
(PC-AT).  The  PC  was  then  used  to  compute  and  suitably 
average  the  power  spectrum  of.vU)  or  the  Poincaie  map 
of  x  and  .v.  For  the  Poincare  map,  the  ADC  was  tin;- 
gered  by  a  pulse  from  the  higher -frequency  one  of  s\u> 
Coinstron  frequency  svnthesi/eis.  The  power  skvtr.i 
were  obtained  from  tieciiiniiiin:  tune  seues.  i  .  tun. 
senes  wtiii  ti  wet,*  ,|l(. 


x=  -  kx  —  x  — J3sin2nx  ! 

+  q ,  sinw,  t  +  q2  sinwt  t  | 

FIG.  2.  Schematic  diagram  of  the  analog  simulator.  The  cir-  i 
cles  marked  (  X )  are  analog  devices  AD534  voltage-multiplier  I 
chips;  the  sine  chip  is  an  AD639;  the  integrators  are  operational  I 
amplifiers  with  feedback  capacitors;  and  summation  is  accom- 1 
plished  by  an  operational  amplifier  with  multiple  input  resistors. 
The  approximately  irrational  forcing  is  accomplished  with  two 
Comstron  frequency  synthesizers  operated  at  /,  =3.5810  kHz 
and  fj  =2.2131  kHz  which  have  a  resolution  of  0.1  Hz  and  a 
short-term  stability  of  a  few  parts  in  10’.  The  time  constant  of 
the  integrators  was  10-4  s.  The  voltages  jc ( r )  and  x(/)  were 
sampled  and  digitized  in  pairs  upon  arrival  of  a  'transistor- 
transistor-Iogic  trigger  pulse  from  the  higher-frequency  syn¬ 
thesizer. 


the  synthesizer. 

It  is  important  that  the  simulator  be  driven  by  signal 
generators  of  higher  resolution  and  better  frequency  sta¬ 
bility  than  is  normally  obtainable  from  ordinary  equip-' 
ment.  Frequency  synthesizers  are  therefore  necessary  in' 
order  that  relative  frequency  drift  not  wash  out  the  spe¬ 
cial  dynamical  effects  being  sought.  Even  so,  the  ratio  ofi 
the  synthesized  frequencies  could  only  approximate  the; 
number  shown  xn  Eq.  (3)  to  the  fourth  decimal  place.  As; 
shown  below,  by  the  agreement  of  our  measurements,' 
with  digital  numerical  simulations,  this  precision  was  evi-i 
dently  sufficient.  The  simulator  was  operated  with  a: 
voltage-scale  factor  of  unity,  which  meant  that  voltage 
measurements  taken  from  it  were  numerically  equivalent, 
to  digital  solutions  of  Eq.  (2).  The  simulator  docs,  how-: 
ever,  scale  time  by  a  factor  equal  to  the  integrator  time' 
constant  r,.  so  that  for  the  digital  simulations,  the  fre¬ 
quencies  must  be  interpreted  as  dimensionless  quantities 
<u'~-«r,.  Since  the  integrator  time  constant  was  always 
10  4  s,  the  dimensionless  frequencies  for  use  in  the  digital! 
simulations  of  Eq.  12)  were  =  2. 2500  and! 
I.  3905  1 0.000 1,  respectively.  It  is  also  to  be  noted  that 
ibe  stability  of  a  few  parts  in  104  is  much  greater  than  the 
fourth  decimal-place  precision  achieved. 

I  lie  sicady-statc  accuracy  of  the  analog  simulator  was 
estimated  to  be  approximately  1%,  and  was  determined 
by  ibe  precision  of  the  resistors  used  Ihc  accuracy  of 
the  A  DM4|.  chip  is  0.25','c.  The  time  scaling  accuracy  of 

r  '  :  was  determined  by  the  precision  of  the  feedback 
•.’•■‘is  trt  the  mtcgialofs.  (nn  simulator  was  operat¬ 


ed  in  a  bandwidth  of  20  kHz,  but  the  largest  frequency  in 
our  data  never  exceeded  5  kHz.  In  order  to  take  advan¬ 
tage  of  the  “'small-signal”  linearity  of  our  operational 
amplifiers,  the  simulator  was  operated  over  an  output  ‘ 
range  of  ±  2  V  instead  of  the  maximum  usable  range  of  ± 

10  V.  In  effect,  the  dynamics  sampled  frequently  only  the 
two  potential  wells  at  si  0.9  V,  as  shown  in  Fig.  1,  and  1 
only  much  less  frequently  visited  the  higher-order  local 
minima.  A  Poincare  map  of  10000  pairs  of  points  could 
be  obtained  in  a  few  seconds.  Speed  is  one  advantage  of 
analog  simulation. 

The  digital  simulations  were  accomplished  using  stan¬ 
dard  techniques.  Specifically,  the  method  used  to  ealeu-  . 
late  the  Lyapunov  spectrum  is  that  due  to  Shimada  and 
Nagashima  and  Benettin  et  a!.  [18].  The  actual  comput¬ 
er  code  used  was  that  due  to  Wolf  et  al.  [19]  with  some 
modifications.  The  Lyapunov  exponent  was  computed 
from  trajectories  and  averaged  until  convergence  to  less 
than  10-4  was  achieved.  For  the  q2—03S  (chaotic)  at-  j 
tractor  3731  trajectories  were  required  for  convergence, 
while  the  ^2::=0.88  (strange  nonchaotic)  attractor  re-  : 
quired  7000.  The  information  dimension  was  obtained  : 
from  the  complete  three-dimensional  phase  space  using  i 
typically  30000  trajectories.  These  numerical  determina-  ! 
tions  of  D i  were  in  good  agreement  with  those  obtained  : 
from  the  two-dimensional  Poincare  sections  in  the  analog 
simulations  with  Djd—  1  =D}d. 
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FIG  '  A  SNCA  lor  A  /!  2  ami  ,h  2.768,  which  aic 

henceforth  identified  as  the  "fixed  conditions”  ami  will)  condo) 
parameter  </.  0  88.  Iiom  the  dynamics  of  Fqs.  (2)  ill  i.e 

Mcasmed  on  the  analog  simulaloi  with  lO(XX)  points  displayed, 
anil  (h)  obtained  hv  iligunl  simulation  lor  ihc  satne  comblioi:-. 
n,,h  10000  poinl>  displaced  I  he  I  v.ipunoc  esptmcnl.  obi  iim  ,1 
■is  dfsctiltvd  in  liir  If  tl,  uas  /  it  |S|)!,U  this  alli.icloi 


A  typical  Poincare  map  of  a  SNCA  for  this  system  is 
shown  in  Fig,  3.  The  attractor  measured  with  the  analog 
simulator  is  shown  in  Fig.  3(a)  and  can  be  compared  to 
one  obtained  for  the  same  conditions  by  digital  simula¬ 
tion  shown  in  Fig.  3(b).  The  attractors  are  quite  similar, 
though  the  digital  simulation  has  fewer  points,  and  the 
fine  scale  of  the  analog  attractor  is  somewhat  obscured  by 
noise.  Nevertheless,  the  agreement  is  quite  convincing. 
Below,  we  discuss  these  and  additional  data  in  more  de¬ 
tail. 

UI.  RESULTS 

Figure  3  shows  a  SNCA  measured  for  4,  =0.88  and 
other  conditions  fixed  as  indicated  in  the  figure  caption. 
Throughout,  we  regard  q2  as  a  control  parameter.  We: 
show,  in  Fig.  4,  a  chaotic  attractor,  again  obtained  from, 
analog  (a)  and  digital  (b)  simulations.  This  attractor  was 
obtained  for  q2~03&,  and  the  digitally  obtained: 
Lyapunov  exponent  was  A=0.156.  At  this  point,  the 
only  criterion  wh’ch  wc  have  used  for  identifying  the  at¬ 
tractor  type,  apart  from  the  obvious  difference  in  appear-: 
ance  (with  the  attractor  in  Fig.  3  appeanng  more  like  the! 
description  [1,3]  “everywhere  discontinuous  but  single: 
valued")  is  the  sign  of  the  Lyapunov  exponent.  We  shall 
now  apply  the  other  criteria  outlined  in  Sec.  I. 

Two  power  spectra,  measured  on  the  analog  simulator 
for  the  fixed  conditions,  are  shown  in  Fig.  5(a)  witn 
q2=0M  and  Fig.  6(a)  with  <72=0.38,  that  is,  for  condi- 
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F  IG.  4.  A  chaotic  attractor  for  the  fixed  conditions  with 
<7.  0.38.  (a)  Measured  on  tha  analog  simulator  1  10*  points), 

and  (b)  by  digital  simulation  (3X10J  points)  from  which 
X  0  1  56  was  obtained. 
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FIG.  5.  (a)  A  power  spectrum  obtained  from  the  analog 
simulator  for  the  fixed  conditions  and  g}  —0.88  (SNCA),  and  <b) 
the  threshold  dependence  of  the  number  of  peaks  N{u)  with 
amplitudes  greater  than  the  threshold  a. 
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FIG  6  (a)  A  power  spectrum  obtained  by  analog  simulation 
for  the  fixed  conditions  and  q,  0  'H  (chaotic),  and  ih)  the 
threshold  dependence  of  the  number  of  peaks 


tion-s  identical  to  those  of  the  attractors  shown  in  Figs  3 
and  4.  The  most  striking  difference  is  that  the  multitude 
of  peaks  evident  in  Fig.  5(a)  are  broadened  in  Fig.  6(a) 
with  only  a  few  dominant  ones  remaining  The  latter 
clearly  results  from  a  chaotic  system,  while  the  former 
might  well  be  supposed  to  result  from  a  very  complex 
quastperiodic  attractor.  We  have  attempted  to  apply  the 
criterion  given  by  Eqs.  (1)  to  these  two  power  spectra, 
that  is,  we  have  counted  the  number  of  peaks  N ( a  s  with 
amplitudes  greater  than  a  threshold  a.  The  results,  plot¬ 
ted  on  a  log-log  scale,  are  shown  below  their  respective 
power  spectra  in  Figs.  5(b)  and  6(b).  While  it  is  not  quite 
possible  to  match  the  results  shown  in  Fig.  5(b)  to  a  pure 
power  law  as  indicated  by  Eq.  1(a),  that  seems  closer  to 
an  appropriate  description  than  that  of  the  results  shown 
in  Fig.  6(b).  It  is  possible  that  the  washing  out  of  the  very 
fine  structure  by  the  inherent  noise  in  the  analog  system 
and/or  the  smearing  of  some  of  the  finer  peaks  in  the  fre¬ 
quency  spectra  by  phase  noise  in  the  synthesizers  may 
have  obscured  the  power-law  dependence  of  N(a )  on  a. 

We  now  turn  our  attention  to  the  information  dimen¬ 
sion,  defined  in  the  usual  way  (6,20): 


D,  —  lim 

f  — o 


/(e) 

log,0(l/e) 


where  e  is  the  box  size  in  phase  space  and  lie)  is  the  in¬ 
formation 

V(e) 

/(<?)=-  2  />, InP,  ,  (6) 

i  -1 

with  Pj  the  frequency  at  which  box  i  is  visited  by  the  dy¬ 
namics.  We  have  measured  this  dimension  on  data  pro¬ 
vided  by  the  analog  simulator  from  the  reduced,  two- 
dimensional  (Poincare)  section.  In  Fig.  7(a)  we  have  plot¬ 
ted  He)  versus  log:o(  l/e)  for  the  two  attractors  shown  in 
Figs.  3  and  4.  Over  the  range  in  e  above  a  lower  limit  set 
by  the  inherent  noise,  we  find  a  good  linear  behavior  of 
lie)  on  log|0(l/<?)  as  shown  by  the  straight  lines,  the 
slopes  of  which  give  D,.  We  find  that  Z)7  =  1. 19  for  the 
attractor  of  Fig.  3  to  be  contrasted  with  D,^  1.58  for  the 
attractor  of  Fig.  4.  The  former  result  is  in  good  agree¬ 
ment  with  both  the  predictions  and  the  numerical  results 
of  Ding,  Grebogi,  and  Ott  [6j  for  a  SNCA.  Since  the 
measurement  of  dimension  on  real  physical  systems, 
which  are  always  to  some  degree  contaminated  by  noise, 
is  often  open  to  criticism,  we  have  tested  our  algorithm 
on  measured  sets  of  data  of  various  sizes  for  the  SNCA  of 
Fig.  3.  These  results  are  shown  in  Fig.  7(b),  where  it  is 
evident  that  data  sets  varying  in  size  from  104  to  105 
points  yield  substantially  the  same  fitted  slope  if  rounding 
at  the  high  and  low  ends  of  the  iog,0(  1  / e )  scale  is  neglect¬ 
ed.  At  the  high  end  there  is  more  rounding  for  the  small¬ 
er  data  sets,  as  expected. 

We  consider  now  the  destruction  of  the  chaotic  attrac¬ 
tor  which  existed  at  q}  0.38  and  its  replacement  by  a 
SNCA  at  a  critical  value  <7?trl,  ‘  0.61 2.  These  results 
were  obtained  by  digital  simulation,  and  the  only  cri¬ 
terion  used  to  determine  the  type  of  dynamics  was  the 
Lyapunov  exponent.  By  measuring  /,  for  each  attractor 
obtained  from  stepping  q,  between  the  two  limits  given 


FIG.  7.  The  information  v 


the  logarithm  of  the  reciprocal 


box  size  obtained  from  analog  measurements  (a)  The  crosses 
are  for  0.88  (SNCA)  and  [he  asterisks  are  for  q:  "0  38 
(chaotic).  The  straight  lines  have  slopes  D,  "'  1.19  (  +  )  and 
£>/Sl.58(,).  The  size  of  the  daia  sets  were  I  O'  points,  (b)  A 
test  of  the  dependence  of  D,  on  die  dala-sel  size  from  analog 
measurements  (  j)  l(f  t  X  )  2  X  10*  (  -  >  4  X  !0\  and  ( .  )  105 
points,  respectively.  The  logarithms  are  base  10. 


above,  we  determined  the  value  of  q\.„  at  which  A  passed 
through  zero.  Treating  Ui  as  an  order  parameter,  we 
then  made  the  plot  shown  in  Fig.  8.  Least-squares  fits  to 
obtain  the  straight  lines  shown  indicate  that  the  results 


FIG  8.  The  1  yapumu  exponent  .is  in  order  parameter  in  the 
destruction  of  chaos  bv  a  SNCA  I  he  results  for  the  chaotic  at 
tractor  are  shown  hv  the  dash-doited  line,  and  tor  the  SNCA  bv 
the  solid  line  Hie  straight  lines  are  least  squares  tils  in  lh"  data 
and  have  slopes  IIM  t  and  0  4M  tor  the  i  liaoln  and  none  hat  it  u 
attractors,  respect  is  civ 


measured  the  two-dimensional  densities  Ptx.x)  using  the 
analog  simulator,  with  the  results  shown  in  Fig  'if a!  for 
the  chaotic  attractor  and  4(b)  for  the  SNCA  We  note 
that  P(x,x)  is  much  more  sharply  defined  for  the  SNCA 
than  for  the  chaotic  attractor.  This  >s  indicated  by  the 
amplitudes  of  the  peaks,  both  densities  having  been  ob¬ 
tained  from  the  same  number  of  digitized  points  (  105). 
The  chaotic  attractor  displays  a  considerably  smaller  am¬ 
plitude  and  shows  less  fine  detail  as  expected  from  a 
higher-dimensional  object.  It  may  be  possible  to  develop 
an  additional  criterion  based  on  the  differences  in  these 
densities. 


FIG.  9.  The  invariant  densities  for  the  attractors  shown  in 
Figs.  3  and  4,  measured  on  the  analog  simulator  using  data  sets  . 
of  10'  points  each:  (a)  for  the  chaotic  attractor  with  the  fixed 
conditions  and  q j  =0.38,  and  (b>  for  the  SNCA  with  fixed  condi¬ 
tions  and  q,  =0.88.  The  velocity  x  is  plotted  on  the  isometric 
axis  vs  the  coordinate  x  on  the  horizontal  axis. 


behave  according  to  a  power  law, 

|'>|0t'  {'HP*  l<?2~?2cn.  I  -  <7> 

with  as0.463  for  the  SNCA  and  asO.  613  for  the  chaot¬ 
ic  attractor.  We  note  also  the  much  greater  variability  of 
A  for  the  SNCA,  an  effect,  noted  in  Ref.  [7],  which  possi¬ 
bly  contributes  to  the  difficulty  encountered  in  attempts 
to  measure  experimentall^/Pfrom  Poincare  sections. 

Finally,  sve  present  also  without  analysis  another  mea¬ 
sure,  the  invariant  density,  which  shows  differences  be¬ 
tween  the  two  attractors  shown  in  Figs.  3  and  4.  We 


IV.  SUMMARY  AND  CONCLUDING  REMARKS 

We  have  demonstrated,  with  measurements  on  an  ana¬ 
log  simulator  —  a  real  physical  system — for  a  particular 
multistable  potential  not  previously  studied,  that  the  dy¬ 
namics  [1-6]  observed  in  bistable  magnetic  beam  experi¬ 
ments  [7,8]  is  robust  to  the  small  inherent  noise  encoun¬ 
tered  in  experimental  systems,  exists  over  a  region  of  pa¬ 
rameter  space  of  nonzero  measure,  and  is  easily  observ¬ 
able  using  criteria  based  on  the  Lyapunov  exponent,  in¬ 
formation  dimension,  and  characteristics  of  the  power 
spectra.  Moreover,  we  have  shown  that  a  SNCA  can  des¬ 
troy  and  replace  a  chaotic  attractor.  We  have  displayed 
this  phenomenon  as  a  classical  phase  transition  using  the 
Lyapunov  exponent  as  an  order  parameter,  whereupon 
we  obtained  substantially  different  power-law  exponents 
for  the  SNCA  compared  to  the  chaotic  attractor.  Final¬ 
ly,  we  suggest  that  our  simulations  can  be  viewed  as  a 
prelude  to  further  experimentation  on  SNCA's  with  irra¬ 
tionally  forced  r.f  SQUID’s,  whose  potential,  with  regard 
to  the  magnetic  flux  internal  to  the  SQUID  loop,  is  simi¬ 
lar  to  the  one  used  here. 
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